THE STOCHASTIC HEAT EQUATION WITH A FRACTIONAL-COLORED 
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Abstract. In this article we consider the stochastic heat equation ut — An = B in (0, T) X R d , with 
vanishing initial conditions, driven by a Gaussian noise B which is fractional in time, with Hurst index 
H 6 (1/2, 1), and colored in space, with spatial covariance given by a function /. Our main result 
gives the necessary and sufficient condition on H for the existence of the process solution. When / 
is the Riesz kernel of order a g (0, d) this condition is H > (d — a) /A, which is a relaxation of the 
condition H > d/A encountered when the noise B is white in space. When / is the Bessel kernel or 
the heat kernel, the condition remains H > d/A. 



1. Introduction and Preliminaries 

Stochastic partial differential equations (s.p.d.e.'s) perturbed by noise terms which bear a "col- 
ored" spatial covariance structure (but remain white in time) have become increasingly popular in the 
recent years, after the fundamental work of [2]. Such an equation can be viewed as a more flexible alter- 
native to a classical s.p.d.e. driven by a space-time white noise, and therefore it can be used to model 
a more complex physical phenomenon which is subject to random perturbations. The major drawback 
of this theory is that it is mathematically more challenging than the classical theory, usually relying 
on techniques from potential analysis. One advantage is that an s.p.d.e. perturbed by a colored noise 
possesses a process solution (under relatively mild conditions on the covariance structure), in contrast 
with its white-noise driven counterpart, for which the solution is well understood only in the sense of 
distributions. Another advantage is the fact that such an equation can lead to a better understanding 
of a complex physical situation. 

This article continues the line of research initiated by [2j , the focus being on a relatively simple 
s.p.d.e., the stochastic heat equation. The novelty comes from the fact that random noise perturbing 
the equation possesses a colored temporal structure given by the covariance of a fractional Brownian 
motion (fBm) , along with the colored spatial structure of [2] . 

We recall that a fBm on the real line is a zero-mean Gaussian process with covariance function 
Unit, s ) = {t 2H + s 2H — \t — s | 2ff )/2, where H £ (0, 1). There is a huge amount of literature dedicated 
to the fBm, due to its mathematical tractability, and its many applications. We refer the reader to [5] 
for a comprehensive review on this subject. 

Recently, the fBM made its entrance in the area of s.p.d.e.'s. We refer, among others, to [5], 
[IB] . [5], [10], [13] or [14]. For example, in the case of the stochastic heat equation driven by a Gaussian 
noise which is fractional in time and with a rather general covariance in space, it has been proved in 
[IB] that if the time variable belongs to [0, T] and the space variable belongs to S 1 (the unit circle) then 
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the solution exists if and only if H > 1/4. The case of the same equation driven by a fractional- white 
noise with space variable in R d has been treated in [8], and it follows that a process solution exists if 
H > d/4. 

We note in passing that very few results are available in the literature, in the case of non-linear 
equations driven by a fractional Gaussian noise, due to difficulties encountered in the stochastic calculus 
associated to this noise. To circumvent these difficulties, a new pathwise method has been developed 
recently in [13] and [14j . treating the nonlinear stochastic heat equation, respectively the nonlinear 
stochastic wave equation. In both these articles, the noise term carries a colored spatial covariance 
structure as well, which is the situation that we investigate in the present paper too. 

In the present article we consider the (linear) stochastic heat equation in the domain [0, T] x R d 
driven by a Gaussian noise B which has a fractional time component, of Hurst index H 6 (1/2, 1), and 
a colored spatial component. Therefore, our work lies at the intersection of the two different lines of 
research mentioned above, namely those developed in [2], respectively |10j . The solution of the equation 
will be given in the mild formulation, but can also be viewed as a distribution solution. Therefore, the 
first step we need to take is to develop a stochastic calculus with respect to the noise B. Since our 
equation is linear, only spaces of deterministic integrands are considered in the present paper. Our 
main result identifies the necessary and sufficient condition on the Hurst index H for the existence of 
the process solution. We should mention that in this case, the solution is a Gaussian process, and hence 
Reproducing Kernel Hilbert Space techniques can be used to investigate its properties. This will be the 
subject of future work. 

In preparation for treating the existence problem in its full generality, we studied first the case of 
the heat equation perturbed by a Gaussian noise which is fractional in time, but white in space (Section 
[2]). In this case, it is known from [8] that the process solution exists if H > d/4, which forces a spatial 
dimension d G {1,2,3}; in the present article, we strengthen this result by proving that the condition 
H > d/4 is in fact necessary for the existence of the solution. In contrast, when a spatial covariance 
structure is embedded in the noise (Section [3]) the condition for the existence of the solution can be 
relaxed so that it imposes no restrictions on the spatial dimension d. When the color in space is given 
by the Riesz kernel of order a, we prove that the necessary and sufficient condition for the existence 
of the solution is H > (d — a)/ 4. This demonstrates that a suitable choice of the spatial covariance 
structure can compensate for the drawbacks of the fractional component. However, it turns out that if 
the spatial covariance is given by the Bessel or the heat kernel, the condition remains H > d/4, whereas 
for the Poisson kernel the condition becomes H > (d+ l)/4. 

This article contains 3 appendices. Appendix A contains a lemma which is heavily used in the 
present paper. This lemma is the tool which allows us to import the Fourier transform techniques from 
R to the bounded domain [0, T]. (So far, this type of techniques have been exploited only on R; see e.g. 
[llj.) Appendix B contains the proof of a technical statement. Appendix C contains a result which is 
essentially due to [2J; we include it since we could not find a direct reference. 

We begin now to introduce the notation that will be used throughout this paper. 

If U C R" is an open set, we denote by T>(U) the space of all infinitely diffcrcntiable functions 
whose support is compact and contained in U. By T>'(U) we denote the set of continuous linear 
functionals on T>(U) which is known as the space of distributions. We let <S(R n ) be the Schwarz space 
of all decreasing functions on R n and <S'(R") be the space of tempered distributions, i.e. continuous 
linear functionals on <S(R n ). For an arbitrary function g on R d the translation by x is denoted by g x , 
i.e. g x {y) — g(x + y). The reflection by zero is denoted by g, i.e. g{x) = g[—x). 
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For any function <fi G <S(M n ) we define its Fourier transform by 

^'PiO — / exp(— i£ ■ x)(f>(x)dx. 

The map T : <S(M") — ► <S(R n ) is an isomorphism which extends uniquely to a unitary isomor- 
phism of ^(M™); this map can also be extended to 6>'(]R n ). We define the convolution (/ * g)(x) — 
J Rn f(x - y)g(y)dy. 

For any interval (a, b) C K and <p £ £2(0, 6), we define the restricted Fourier transform of y> 
with respect to (a, 6) by: 




ip{t)dt. 



As in [12] , if / G £i(0, T) and a > 0, we define the fractional integral of / of order a by: 

= ^ f \u~tT- l f{u)du. 
Finally, we denote by Bt,(M. d ) the class of all bounded Borel sets in ~R d . 

2. The Fractional- White Noise 

The purpose of this section is to identify the necessary and sufficient condition for the existence 
of a distribution-solution of the stochastic heat equation, driven by a Gaussian noise which is fractional 
in time and white in space. 

Our main result is similar to Theorem 11, [2], which identifies the necessary and sufficient 
condition for the existence of a distribution-solution of an arbitrary s.p.d.e.'s driven by a Gaussian 
noise which is white in time and colored in space. 

In the first subsection we examine some spaces of deterministic integrands, which are relevant 
for the stochastic calculus with respect to fractional processes, and we explore the connection with 
Dalang's theory of stochastic integration. In the second subsection, we describe the Gaussian noise and 
its stochastic integral. In the third subsection, we introduce the process-solution and the distribution- 
solution of the stochastic heat equation driven by this noise, and we identify the necessary and sufficient 
condition for the existence of these processes. 

2.1. Spaces of deterministic integrands. We begin by introducing the usual spaces associated 
with the fractional temporal noise. Throughout this article we suppose that H £ (1/2,1) and we let 
a H = H(2H - 1). 

Let 7i(0, T) be the completion of T>(0, T) with respect to the inner product 

( ( P,i J )H(o,T) = ct H / (p{u)\u~v\ 2H 2 ip{v)dvdu 
Jo Jo 

= a H c H [ fo.T^rJV^WIrr^-^r, 
Js. 

where the second equality follows by Lemma [All (b) with c H = [2 2 ^- If >ir 1 / 2 ]- 1 r(H - l/2)/T(l - H). 
Note that cjj = cfoH-i, where the constant q a is defined in Lemma lA.ll (Appendix A). 

Let £(Q,T) be the class of all linear combinations of indicator functions l[o,t]>i G [0,T]. One 
can see that 7i(0,T) is also the completion of £(0,T) with respect to the inner product 



(l[o,t] ; 1[o,«])w(o,t) — Rn{t, s). 
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(To see this, note that every l[o, a ] € £(0, T) there exists a sequence ((p n )n C T>(0, T) such that <p n (t) — > 
l[o.a](^))Vt and supp C if for all n, where if C (0, T) is a compact. By the dominated convergence 
theorem, it follows that \\(p n — ho,t] II h(o,t) ~ * 0.) 

In is important to emphasize that the space H(0,T) may contain distributions. We will jus- 
tify this statement, using the argument of [2]. First, note that Tt(0,T) c H(R), where H(R) is the 
completion of 2?(R) with respect to the inner product 

2H-2, 



(v> V')'H(iR) = I I ip(u)ip(v)\u - v\ dvdu. 

JR JR 

The space 7i(R) appears in several papers treating colored noises. (In fact W(R) is a particular instance 
of the space v\fl of [I], in the case /i(d£) = |£|-( 2H - x )d£ and d = 1.) From p. 9 of [2], we know that 

H(R) C W(R) := {5 € S'(R);.FS is a function, f \TS{t)\ 2 \t\- (2H ^ 1] dr < oo}. (1) 

Jr 

Since |r| 2 < 1 + |t| 2 , one can easily see that 7?(R) C U^ {H ^ 1/2) (R), where 

w -(H-i/ 2 ) (M) . = | 5 e 5'(R) ; ^5 is a function, / |^5(r)| 2 (l + | ^ 1 2 ^ - (J^- 1/2) < ^ 

is the fractional Sobolev space of index —(H~ 1/2) (seep. 191, [6]). Therefore, the elements of H(R) are 
tempered distributions on R of negative order —{H — 1/2). (This was also noticed by several authors; 
see e.g. p. 9, [5], or [IT], [H]) 

On the other hand, similarly to (jTJ) , one can show that 

H(0,T) C H(0,T) := {S e S'(R); .Fo.tS is a function, / |.Fo,TS(T)|Vr (2H_1) dT < 00}, (2) 

where £/ie restricted Fourier transform Tq^tS of a distribution S G <S'(R) is defined by: (Tq^tS,^) = 



Remark 2.1. One can show that (see p. 10, [5]) 

L 2 (0,T) c Li/ir(0,r) c \H(0,T)\ c W(0,T), (3) 

where \H(0,T)\ = {/ : [0, T] xl^I measurable; J T f Q T |/(u)| \f(v)\ \ \u - v\ 2H - 2 dudv < 00}. 

A different approach of characterizing the space 7i(0, T) is based on the transfer operator. We 
recall that the kernel Kn(t, s),t > s of the fractional covariance function Rh is defined by: 

rt 



K H (t, s) = c^ 1 / 2 ^ y (u - a ) ff - 3 /a u »-V2 dU) 



where c* H = {a H T(3/2 - H)/[T(2 - 2H)T(H - 1/2)]} 1 / 2 . 

Note that Rn(t,s) = J* As KH(t,u)K(s,u)du (see p. 7-8, [9]) and hence 



(^lf 1 [o,t],if|fl[o, s ])L2(o.T) = / K H (t,u)K H (s,u)du = R H (t,s) = (l[o,t], 1[o,s])«(o,t), 

Jo 

i.e. if|f is an isometry between (f(0,T),(',')^( ,T)) and £2(0, T). Since H(0,T) is the completion of 
£(0,T) with respect to (•, ■)'H(o,t)i this isometry can be extended to 7i(0,T). We denote this extension 
by K* H . In fact, one can prove that the map K* H : 7i(0, T) — > £2(0, T) is also surjective (see the proof 
of Lemma [~~ 
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We are now introducing the space of deterministic integrands associated with the a noise which 
is fractional in time and white in space. This space was also considered in [lOj . 

More precisely, let TL be the completion of 2?((0, T) x R d ) with respect to the inner product 

{<p,tp)u = &h / / / p(u,x)\u — v\ 2H ~ 2 ip(v,x)dxdvdu 
Jo Jo Js. d 

= a H CH / / J r o.Tip(T,x)J 7 o,Tip(T,x)\T\^ (2H ^ 1 ' > dxdT 

Jr JR d 

JR d 

where the second equality above follows by Lemma lA. 11 (b). and the third is due to Fubini's theorem. 

If we let £ be the space of all linear combinations of indicator functions l[o,t]xA> t G [0, T], A G 
Bb(M. d ), then one can prove that TL is also the completion of £ with respect to the inner product 

(l[o,t]xA, 1[o,s]xb)h = Rn(t, s)X(A n B) 

where A is the Lebesgue measure on R d . (The argument is similar to the one used in the temporal 
case.) Similarly to (J2J), one can show that 

TL C H := {S : R d -> S'(R); F ,tS(-, x) is a function Vie G R d , (t, x) i-> Tq,tS{t, x) is measurable, 

and f f \F .TS{T,x)\ 2 \T\- {2H - 1] dxdT < oo}. 

Jr jR d 

Using ((3]) and the fact that 

||5||^= / \\S(;x)\\ 2 n{0T) dx, VSeW, 

jR d 

one can show that 

L a ((0, T) x R d ) c\TL\(lTL(Z L 2 (R d ; TL{0, T)), 
where |7i| = {</?: [0, T] x R d measurable; J T J T J Rd \<p(u, x)\\ip(v, x)\\u — v\ 2H ~ 2 dxdvdu < oo}. 

The next result gives an alternative criterion for verifying that a function ip lies in TL. (It can 
be compared to Theorems 2 and 3 of [2].) 

Theorem 2.2. Let <p : [0,T] x M. d — > R be a function which satisfies the following conditions: 

(i) p{-,x) G L 2 (0, T) /or ewn/ a; G R rf ; 

(ii) (t, x) h- > To,tp{t,x) is measurable; 

( iU ) 4Jk<< I-^o.t^Ct,^)! 2 ^!-^- 1 )^^ < oo. 
TTien tp TL. 

Proof: Similarly to Proposition 3.3, [IT], we let 

A = {ip : [0, T] x R rf -> R; or) G L 2 (0, T) Vx, (r, x) Fq,t<p(t, x) is measurable, and 

IMl! :=c i/ / l^o,T^(r, : r)| 2 |T|-( 2H - 1 )cfa ; dr<oo} 
Jr jR d 

A = {p:[0,T}xR d ^R; \\p\\ 2 A :=c 2 [ [ [l"l 1/2 (u H - 1/2 p(u, x))(s)} 2 s^ 211 ^ dxds < oo} 

where ci = ancn and c 2 = {cflT(if — 1/2)} 2 . We now prove that 

Ac A and = |M|a, V<y9 G A. (4) 
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Let ip £ A be arbitrary. Since (fi{-,x) £ L 2 (0, T) C H(0, T) and K* H is an isometry from H(0, T) 
to L 2 (0, T), we have \\<p(-, x)\\^ Q T j = \\K^jtp(-, £)||| 2 ( T \ for all x £ R d , that is 

ci / |^o^(T,x)| 2 |r|-( 2H - 1 )dr = C2 /" T [7^r 1/a (u ir - 1/ V(«,a:))(*)] 2 a- (2H - 1) d«, Vx e R d . 



Integrating with respect to dx and using Fubini's theorem, we get ||</?|| A = ||</?||a- The fact that 
||</j|| a < oo forces ||v||a < oo, i.e. 93 £ A. This concludes the proof of (H}. 
Next we prove that 

£ is dense in A with respect to || • ||a- (5) 

Let ip £ A and e > be arbitrary. Since the map (s,x) 1— ► 1 ' 2 (u H ~ 1 / 2 (p(u, x))(s) belongs 
to Z,2((0, T) x R d ,(iA_f/ x c?x) where Xh(s) = s~( 2H ~ l 'ds, there exists a simple function g(s,x) = 
ELi ^Ih.rftjfsJUiW on (0,T) x K d , with (< t el 1 0<c l <4<r and A k C M d Borel set, such 
that 

/ T / [7^T 1/2 (u /f - 1 /VKa;))(s)-5(s ! a;)] 2 s-^- 1 )da;ds<£. (6) 

JO JT$. d 

By relation (8.1) of [12], we know that there exists an elementary function l k £ £(0,T) such that 



[W)M - i^: l/2 {u H -^i k {u)){s)f s -^ H -^d s < e/c g , 

where we chose C g := nJ2k=i tfc^iAk)- Define l(s,x) = bklk(t)lA k (x) and note that I £ £. Then 

[g{s, x) - I*Z Xl2 {u H - x l 2 l{ Ul x)){s)fs-^ H -^ds < e. (7) 



From © and ((7]), we get 

[ T f [I^: 1/2 (u H -^ 2 ^(u,x))(s) -I^i^-^l^x^fs-^-^dxds < 4e, 
Jo Js. d 

i.e. \\ip — l\\\ < 4ec2- This concludes the proof of . 

From ([4]) and ([5]), we infer immediately that £ is dense in A with respect to || ■ || A . Since 
II ' II A = II ' lift an d H is the completion of £ with respect to || • it follows that A C TL. This concludes 
the proof of the theorem. M 

As in jTD], we define the transfer operator by: 

( K H 1 [o,t]xA)(s,x) := K H (t, s)l[ 0ii ] Xjl (s,x). (8) 

Note that 

( K H^[0.t]xA, ^lfl[0, ;S ]xB)L 2 ((0,T)xR' i ) = (J K H (t, u)K H (s , u)dl^J (1a, Lb) L 2 (R d ) 

= R H (t, s)X(A n B) = (l[ .t]xA, l[o,s]xfl)«i 

i.e. -ftfjjj is an isometry between (£, (•, -)^) and L2((0,T) x R d ). Since 7i is the completion of £ with 
respect to (■,•)«, this isometry can be extended to TC. We denote this extension by K^. 

Lemma 2.3. : H — > L 2 ((0,T) x R d ) is surjective. 
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Proof: Note that l [0 , t]xA £ K* H {H) for all t £ [0, T],A£ B b (R d ); hence ip £ K^(H) for every ip £ £. 
Let / £ L2((0, T) x M d ) be arbitrary. Since £ is dense in ^((0, T) x M. d ), there exists a sequence 
(fn) n C £ such that || /„ - /||L 2 ((o,T)xR d ) — * 0. Since /„ £ Kft(H), there exists tp„ £ H such that 
/„ = Kftip n . The sequence ((p n )n is Cauchy in Ti: 

\\fn - <Pm\\n = W^H^n ~ K^Prn | L 2 ((0,T) xM d ) = \\fn ~ fm I L 2 ((0,T) xM d ) ~ ► 

as m,n ^ oo. Since TL is complete, there exists ip £ 7i such that \\<p n — p\\u — > 0. Hence ||/„ — 
^H¥>lk 2 ((o,T)xR-*) = II^H^n-^-H < /?llL 2 ((o,T)xR<') -> 0. But || /„ - /|| £ 2 ( (0 ,t) xm") -> °- We conclude that 
= /• ■ 



Remark 2.4. Note that for every p> £ £ 

/ T FlKv C T /r\ H - 1 / 2 

f(r,x)^-(r, S )dr = c* H J <p(r,x){~) {r- s ) H ^l 2 dr 

= 4r(ff-^r(^)j^(/-vv(^))( s ). (9) 

Using Lemma 12.31 we can formally say that 

H = {ip such that (s,x) i ► s-( H - 1 / 2 )/^T 1/2 (u H - 1 / 2 ¥ .(u,a;))(s) is in L 3 ((0,T) x 

2.2. The Noise and the Stochastic Integral. In this paragraph we describe the Gaussian noise 
which is randomly perturbing the heat equation. This noise is assumed to be fractional in time and 
white in space and was also considered by other authors (see [TO]). 

Let F = {F(ip); ip £ 2?((0, T) x R d )} be a zero-mean Gaussian process with covariance 

E(F(<p)F(tP)) = {<p,ip) n - (10) 

Let H F be the Gaussian space of F, i.e. the closed linear span of {F(ip); ip £ V((0, T) x in L 2 (fl). 

For every indicator function Im^xA S £, there exists a sequence (<p n ) n C T>((0,T) x R d ) 
such that ip n — > l[o,t]xA an d supp ip n C K, Vn, where K C (0, T) x R d is a compact. Hence \\p n — 
l[o,i]xA||w — *■ and E(F((p m ) - F(tp n )) 2 = \\p m ~<Pn\\H — > as m, n — > oo, i.e. the sequence {F(v?„)}„ 
is Cauchy in Z^^)- A standard argument shows that its limit does not depend on {ip n } n - We set 
F t (A) — F(l[o,(] x j) =L 2 (Q) hm„F((p„) 6 H F . We extend F by linearity to £. A limiting argument 
and relation (fTD} shows that 

E(F(<p)F(ip)) = {<p,Tp)n, V<P,*l>e£, 

i.e. y> i— > F(y) is an isometry between (£, (•, and H F . Since Tt is the completion of £ with respect 
to (-, -)n, this isometry can be extended to H, giving us the stochastic integral with respect to F. We 
will use the notation 

F(<p) = f f <p(t,x)F(dt,dx). 

JO JR d 

Remark 2.5. One can use the transfer operator to explore the relationship between F(ip) and 
Walsh's stochastic integral (introduced in [17] ). More precisely, using Lemma [2~3l we define 

W{4>) :=F((^)- 1 ( ( /»)), cj ) £L 2 ((0,T)xR d ). (11) 

Note that 

E(W(<t>)W(ti)) = ((^)- 1 (^),(^)- 1 (t ? )) w = (4>,v)l 2 ((o,t)^), 
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i.e. W — {W(<fi);<fi G ^((0, T) x M. d )} is a space-time white noise. Using the stochastic integral 
notation, we write W(<p) = ^ J Rd (f>(t,x)W(dt,dx) for all G L 2 ({0,T) x R d ). (Note that W(<f>) is 
Walsh's stochastic integral with respect to the noise W.) 

Let H w be the Gaussian space of W, i.e. the closed linear span of {W((j>); 4> G ^((0, T) x R d )} 
in L<2,(Vt). By (jTT|) , we can see that H w — H F . The following diagram summarizes these facts: 



K. 



n 



H 



F 



H 1 




L 2 ((0,T)xK rf ) 

F{fp) = W{K* hV ), Vtp G H, i.e. 

Jo /r- = Jo T / H -(^wV)(*. ^^(dt, dx), G W. 



In particular, F(t, A) — J Q f A Kf{(t, s)W(ds,dy). This relationship will not be used in the present 
paper. 

2.3. The Solution of the Stochastic Heat Equation. We consider the stochastic heat equation 
driven by the noise F, written formally as: 



v t — Av = F, in (0, T) x R d , v(0, •) = 0, 

where Aw denotes the Laplacian of v, and v t is the partial derivative with respect to t. 
Let G be the fundamental solution of the classical heat equation, i.e. 



G(t,x) 



(4vrf)- d / 2 cxp (-^rj ift>0,xGi? d 
if t < 0, x G R d 



(12) 



(13) 



Let g tx (s,y) :— (G tx ) ( s ,y) — G(t — s,x — y). The following result shows that the kernel G has the 
desired regularity, which is needed to apply Theorem 12.21 

Lemma 2.6. Iftp = r)*G, where 77 G X>(0, T) x W l ), then 

ip(-,x) Gi 2 (o,T) \fxem d . 



Proof: Without loss of generality, we suppose that f](t,x) = (f>{t)ip(x), where <f> G T>(0,T) and ip G 
"\ Using Minkowski's inequality for integrals (see p. 271, [15]), we have 



1/2 



\ip(t, x)\ 2 dt 



< 




(s)ip(y)G(s -t,y- x)l {s>t} dyds 



1/2 



dt 



\<i>{s)i>{y)\ ( l o \G{s-t,y-x)\ z dt) d,,ds. 



1/2 



Using the change of variables s — t = t' and 1/t' = u, we get 

1 



|G(s-t,y- :r)| 2 di 



00 r 1 

d-3 e -|,,-*|V2 du < 
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and 

/ _ \ 1/2 rp 

I \v(t,x)\ 2 dt) <C [ \<j>(s)\ [ ^ y }}_ 1 dyds<C(\i>\*R 1 )(x)<o , 

o j Jo jR d \v ~ x \ 

where Ri(x) — 71, dM - ^ -1 -* is the Riesz kernel of order 1 in K d and the convolution \ip\ * R\ is well- 
defined by Theorem V.l, p. 119, [HJ. ■ 

The next result gives the necessary and sufficient condition for the existence of the process 

solution. 

Theorem 2.7. If 

H>±, (14) 

then: (a) g tx € \H\ for every (t, x) G [0, T] x W l ; (b) r/*G eH for every rj e P((0, T) x R d ). Moreover, 
\\gtx\\H < oo\/{t,x) e [0,T}xR d if and only if {Tft holds. (15) 
Remark 2.8. Note that condition (I14|) implies that d € {1,2,3}, since H < 1. 

Proof: (a) We will apply Theorem 12.21 to the function g tx . Note that for every fixed y 6 K d 

fT ft -, foo 

\g tx (s,y)\ 2 ds = C / 7 _L_ e -l»-*l a /[2(*-.)] ds = C / u d -\-\y-^du < oo, 
i-e. gtx{-,y) S 1/2(0, T). Clearly (r, y) i— > J r a l Tgtx(' r , y) is measurable. We now calculate 



For this, we write 

htxfn = a H c H I \r\ - ,:i "-- n / (/ r —v, ,.(,s.//)r/.w I f / < IT (/'. i/Ulr I ili/tlr 



i-T \ I r T 

|-(2i?-l) 



a H c H 

ft r t 



t r t 



^0 



| r |-(2iT-l) / / e -*r(.-r)j( rj s ) drdsdT , 



= a H f f \s-r\ 2H - 2 I(r,s)drds 
Jo Jo 

where we used Lemma A.l (Appendix A) for the last equality and we denoted 

I(r, s) := { g tx (s, y)g tx (r, y)dy = (2t - s - r)^ 2 



We obtain that 



\\gtx\\n = a H 



f f \s-r\ 2H - 2 (2t-s-r)- d / 2 drds = a H ( f \u - v\ 2R - 2 {u + v)- d / 2 drds. 
o Jo Jo Jo 



Relation (fT5|) follows, since the last integral is finite if and only if 2H > d/2. In this case, we 
have = \\9tx\\n < 00 ( since 9tx > 0), and hence g tx £ \H\. 

(b) We will apply Theorem 12.21 to the function <p = rj * G, since if(-,x) E £2(0, T) for every 
x S R d , by Lemma \2. 61 By writing 

Fo,t<p{t,x) = / e~ lTS / T)(u + s,y)G(u,y - x)dyduds, (16) 



o J«, d 
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we see that (r, x) t— > Fo.TPir, x) is measurable. We now calculate 

\\<p\\ 2 H :=a H c H f f \F 0<T <p.(T,x)\ 2 \T\-( 2H -VdxdT. 

Jr JR d 

Using (fl6|) , we get 

IMItt = a B c B f H"^" 1 ) f f T e— ^ f f [ T * [ T 7 ' r,(u + 8,y)ri(y + r,z) 

JR JO JO JR d JR d Jo Jo 

J(u, v, y, z)dv du dz dy dr ds dr 

pT />T /• /> />T—s />T—r 

= a H / / \s-r\ 2H ~ 2 / / / / r/(u + s,y)n(v + r,z) 

Jo Jo JR d JR d Jo Jo 

J{u,v,y, z)dv dudz dy dr ds, (17) 
where we used Lemma IA.1I (Appendix A) for the second equality and we denoted 

J(u,v,y,z) := [ G(u,y-x)G{v,z-x)dx = exp\-^- — (u + vy d/2 . 
jR d { 4(u + V) J 

Clearly 

J(u,v,y,z) < {u + i>y d/2 (18) 
for every u G (0, T - s),v G (0, T - r). Using dT7J), (HHJ) and the fact that n G D((0, T) x R d ), we get 

r-T pT p p pT—s pT — r 

\W\\h < oc H \s~r\ 2H - 2 / / / \r)(u + s,y)r]{v + r,z)\ 

Jo Jo JR d JR d Jo Jo 

(u + v)~ d ^ 2 dv du dz dy dr ds 

nT nT r T ~ s r T - r 

< a H C v / / \s-r\ 2H - 2 / / {u + v)- d/2 dvdudrds 
Jo Jo Jo Jo 

T r T 



JO 



< a H C v / / (u + v)- d / 2 {T-u) 2H dvdu, 



where for the last inequality we used Fubini's theorem and the fact that f Q " f Q v \s — r\ 2H 2 dr ds = 
R H (T-u,T-v) = [{T-u) 2H + {T-v) 2H -{u-v) 2H ]/2. The last integral is clearly finite if 2H > d/2, 
i.e. H > d/4. M 

Under the conditions of Theorem 12. 7i F(g tx ) and F(r) * G) are well-defined for every (t,x), 
respectively for every r\ G 2?((0, T) x Mr) and we can introduce the following definition: 

Definition 2.9. a) The process {v(t,x);t € [0,T],x G R d } defined by 

v(t,x) :=F(g tx )= [ [ G(t - s,x - y)F(ds,dy) (19) 



is called the process solution of the stochastic heat equation (I12| . with vanishing initial conditions, 
b) The process {v(r));r] G V((0,T) x R d )} defined by 

v{m):=F(r)*G)= ( [If [ r){t + s,x + y)G(s,y)dyds) F(dt,dx) 



is called the distribution- valued solution of the stochastic heat equation (jT2]l . with vanishing initial 
conditions. 

Lemma 2.10. The process {v(t,x);t G [Q,T],x G R d } is L 2 (Q)- continuous. 
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Proof: We first prove the continuity in t. We have 

rt+h 



E\v(t + h,x)- v(t,x)\ 2 = E 



9t+h,x(s,y)F(ds,dy) 



JTBL d 



9tx(s,y)F(ds,dy) 



E 



< 2E 



JR d 
t 



(gt+h, x (s, y) - gtx(s, y))F{ds, dy) + 



(gt+h,x(s, y) - gtx(s, y))F(ds, dy) 



t+h 



2E 



gt+h,x(s,y)F(ds,dy) 



t+h 



gt+h,x(s,y)F(ds,dy) 



:= 2I 1 (h) + 2I 2 (h). 
We first treat the term I\(h). Note that 

h(h) = \\F((g t+h ,x -5ta:)l[o,t])||£ 2 (n) = \\(9t+h,x - 9tx)\o,t}\\ 



H 



o Jo 



(gt+h,x - gtx)(u, y)\u - v\ (gt+h,x - 9tx)(v, y)dy dv du. 



The continuity of G(t,x) with respect to t shows that the integrand converges to zero as h — ► 0. By 
applying the dominated convergence theorem (using the fact that HsfeH'H < °°), we conclude that 
h(h) -> 0. 

For the term 12(h), we have 

h(h) = \\F{gt+h,x^[t,t+h])\\L 2 (ci) = \\9t+h^[t,t+h]\\n = 

r t+h r t+h r 

gt+h,x(u,y)g t +h,x(v,y)\u - v\ 2H ~ 2 dydvdu 



OLH 



t-h J t-h 



gtx(u',y)gtx(v',y)\u' - v'\ 2 2 dydv' du 



Since l(t-h,t)^(t-h,t) ~ * as h — » and H'fell'H < °°: we conclude that 12(h) — ► 0, by the dominated 
convergence theorem. 

We now prove the continuity in x. We have 



E\v(t,x + h) - v(t,x)\ 2 = E 



(gt,x+h(s, y) - g tx (s, y))F(ds, dy) 



\(gt,x+h - gtx)\o,t]\\n 



= a H 



Jo 



(gt,x+h - 9tx)(u,y)\\u - v\ (gt,x+h - g tx )(v,y)dydvdu. 



By the continuity in x of the function G(t, x), the integrand converges to as h — > 0. By the dominated 
convergence theorem, we conclude that E\v(t, x + h) — v(t, x)\ 2 — ► as h — ► 0. ■ 



Theorem 2.11. Let {v(r));i] S T)((0,T) x M. d )} be the distribution-valued solution of the stochastic heat 
equation \V$l . 

In order that there exists a jointly measurable and locally mean-square bounded process Y = 
{Y(t,x);t£ [0,T],x€ K rf } such that 

fT 



Ufa) 



Y(t,x)n(t,x)dxdt V77 E T>((0, T) x R d ) a., 



(20) 



it is necessary and sufficient that holds. In this case, Y is a modification of the process v 

{v(t,x);t£ [0,T],XE R d } defined by £W\). 
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Proof: The necessity part is similar to the proof of Theorem 11, [2] . Suppose that there exists a jointly 
measurable and locally mean-square bounded process Y = {Y(t,x);t E [0, T],x <G M. d } such that (f2"0|) 
holds. 

Let (to, %o) S [0, T] x R d be fixed. Let A, ip be nonnegative functions such that A € T>(0, T),ip G 
2?(R d ) and \(t)dt — J Rd tp(x)dx = 1. Set A„(t) = n\(nt) and ip n (x) = n d ip(nx). Define r) n (t,x) = 
K(t - t Q )ifj n (x - x ). 

We calculate iS|i;(77 n )| 2 in two ways. First, using (|20|) . we get 

-E|w(?7n)| 2 = / E\Y(t,x)Y(s,y)\ rj n (t,x)ri n (s,y)dydxdsdt. 

J[0,T]xR d 

Using Lebesgue differentiation theorem (see Exercise 2, Chapter 7, [H]), we get 

\imE\v{ Vn )\ 2 = E\Y(t ,x Q )\ 2 . (21) 

Secondly, 

E\v( Vn )\ 2 = E\F( Vn *G)\ 2 = \\vn*G\\ 2 n = I f \T ,T(Vn*G)(r,x)\ 2 \T\- ( - 2H - 1 Urdx. 

Js. d Jm 

We claim that, for every r£l,ie M. d , we have: (see Appendix B for the proof) 

limJ r 0j T(??n * G)(t,x) = To,T9t x a {T,x). (22) 

n 

Using Fatou's lemma, ([21]) and ([22]) . we get 

= / / \^.T9t 0X0 (r,x)\ 2 \r\^ 2H -^dTdx 

< Ik/ ( \T Q , T {rin*G){T,x)\ 2 \T\- {2H -^dTdx 
= IjrnEH^ 2 = E\Y(t 0l xu)\ 2 < oo, 



l5t a;ollw 



which forces iJ > rf/4, by virtue of (fT5l) . 

We now prove the sufficiency part. By Lemma [2.101 we know that the process v defined by 
(Tr^)) is continuous in L 2 (£l). Hence, it is continuous in probability and by Theorem IV. 30, [I], there 
exists a jointly measurable modification Y of v. We have 

E(Y(t, x)Y(s, y)) = E(v(t, x)v(s, y)) = E(F(g tx )F(g sy )) = (g tx , g sy ) H (23) 

and 

E(v(r,)Y(t, x)) = E(v(r,)v(t, x)) = i?(F(ry * G)F(<? 4x )) = (rj * G, g tx ) n . (24) 
To prove that (|2"0|) holds, we will show that 

2 / 



/ Y(t,x)r)(t,x)dxdt 

o il 1 * 
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By Fubini's theorem and ([23]) . we get 

T 2 

"J p 

E I I Y(t,x)r)(t,x)dxdt 

JR d 



r](t, x)r)(s, y) (g tx ,g S y)Hdy dx ds dt 

([0,T]xR d ) 2 



7](t,x)T](s,y) ( / / To,T9tx{T, z)J : 0tT g sy (T 7 z) \t\ {2H ^ dzdr) dydxdsdt 

([O.T]xR d ) 2 \JRJR d / 

\F ,t(v * G)(t, z)\ 2 \T\~ (2R - 1] dzdT = \\n* G\\ 2 H = E\F( V * G)\ 2 = E\v(r,)\ 2 . 



On the other hand, using Fubini's theorem and (|24| . we get 

E\v(ri) / / Y(t, x)r)(t, x)dxdt \ = / r)(t,x){r] * G, g tx )ndxdt = 

y Jo y i[o,T]x« j 

/ r?(t, x) ( f [ T ,t{v * G)(t, z)F ,Tgtx(T,z) \t\~ {2H - 1] dzdT ) rixdi = 

JlO.TlxR^ VJRJR 1 ' / 

To, T (n * G)(r, z)! 2 ^!-^- 1 )^ - E\v( V )\ 2 . 



3. The Fractional-Colored Noise 

In this section we examine the process-solution and the distribution-solution of the stochastic 
heat equation driven by a Gaussian noise which is fractional in time and colored in space. Most of 
the results of this section are obtained by mixing some colored spatial techniques with the fractional 
temporal techniques of Section [2] The results of this section can therefore be viewed as generalizations 
of the results of Section [2j The details are highly non-trivial. 

The structure of this section is similar to that of Section [21 We first describe the spaces of 
deterministic integrands, then we introduce the Gaussian noise and the associated stochastic integral, 
and finally we examine the solution of the stochastic heat equation driven by this type of noise. 

3.1. Spaces of Deterministic Integrands. We begin by introducing the space of deterministic inte- 
grands on M. d . 

We say that a function / : R — > K is a spatial covariance function, if it is the Fourier transform 
of a tempered measure [i on R d , i.e. f(x) — J Rd e^ l ^' x ^i(d£ l ). 

Let V(M. d ) be the completion of T>(M. d ) with respect to the inner product 



i t P,ip)v(R d ) = / f(x)f(x-y)tp(y)dydx= .^(O^KMCh 



where •^vCC) : — Jr<i e~ l ^ x ip(x)dx denotes the Fourier transform with respect to the x-variable. 

Equivalently, we can say that V(M. d ) is the completion of £ (K d ) with respect to (•, -)-p(R d )i where 
£(M. d ) is the space of all linear combinations of indicator functions 1a(x), A S Bb(M. d )) 

The basic example of spatial covariance function is f(x) = S(x), which gives rise to the spatial 
white noise. More interesting covariance structures are provided by potential analysis. Here arc some 
examples (see e.g. p. 149-151, [B], or p. 117-132, [T5]; our constants are slightly different than those given 
in these references, since our definition of the Fourier transform does not have the 2tt factor): 



14 



RALUCA M. BALAN AND CIPRIAN A. TUDOR 



Example 3.1. The Riesz kernel of order a: 

f(x) = R a (x):=j a , d \x\- d+a , 0<a<d, 
where la4 = 2 d - a Tt d l 2 T{{d - a)/2)/T{a/2). In this case, n(d£) = \C\~ a d^. 
Example 3.2. The Bessel kernel of order a: 

/•OO 

f(x) = B a (x) := i a / ^"-^-le-e-HVMrf^ Q > , 



where i a = (47r)"/ 2 r(a/2). In this case, n(d£) = (1 + |£| 2 )~ q/2 A; and V(R d ) coincides with H~ a/2 (R d ), 
the fractional Sobolev space of order —a/2; see e.g. p. 191, [6|. 

Example 3.3. The heat kernel 

f{x)^G a {x)^ 1 l d e-\^'^\ a>0, 
where d = (47ra)~ d/2 . In this case, /z(c/£) = e^ 2 "^ 2 ^. 
Example 3.4. The Poisson kernel 

f(x) = P a (x) := 7 - d (|x| 2 + a 2 )"( d+1 )/ 2 , a > 0, 
where = 7r-( d+1 )/ 2 r((d + l)/2)a. In this case, fi(d£) = e^^d^. 

Remark 3.5. The space V defined as the completion of 2?((0, T) x M. d ) (or the completion of £ ) with 
respect to the inner product 

(ip,ip) v = / / tp(t,x)f(x-y)ip(t,y)dydxdt= (<p(V)^(V))p(i")*- 
Jo Jm d JR d Jo 

has been studied by several authors in connection with a Gausssian noise which is white in time and 
colored in space. One can prove that V C L 2 ((0, T); 7 3 (M d )); see e.g. [2], or pp. 

In what follows, we need to extend the definition of V(M. d ) to allow for complex- valued functions. 
More precisely, let T>c(R d ) be the space of all infinitely differentiable functions ip : M. d — > C with compact 
support, and T'c{R d ) be the completion of T>c{R d ) with respect to 

(^»-p c (R d ) = / / <fi(x)f{x -y)tp(y)dydx. 



Since V(R d ) C T>c(R d ) and (tp,ip)'p(M. d ) = {Vi^VciX*) ^ or evei T <= D(R d ), we conclude that 

V{R d ) C V c (R d ). 

We are now introducing the space of deterministic integrands associated with a Gaussian noise 
which is fractional in time and colored in space. This space seems to be new in the literature. More 
precisely, let TtV be the completion of T>((Q,T) x R d ) with respect to the inner product 



{<p,ip)nv — a H I / / / <p{u,x)\u - v\ 2H 2 f(x-y)ip(u,y)dydxdvdu 

Jo Jo JR d J~BL d 

= a H c H / \t\~@ h ~ 1 ') / / f(x - y)J 7 o_ T ip(T,x)J :: o,Tip(T,y)dydxdT 

Jr JR d JR d 



anc-H \ \r 



-(2H 1) (jr 0jT( ^( T5 .),F T tP(t, ,)) v ( R d)dr, 
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where the second equality follows by Lemma I A. II (Appendix A). In particular, 

IMIwp = a HCH / ||^0,Tp(T,Ollp c (tt<«)M _{2H_1 W, 



or equivalently |M|^> = / R \\T ^{t, OII^R^Aff (dr), where Air(dr) = a ff c H |T| < - 2H ^dr. 

One can prove that 7iP is also the completion of E with respect to the inner product 

(l[0,t]xAi 1[0,s]xb)hP = #J?(*, s)(1a, ls)-p(Rd)- 

Clearly |WP| C WP, where \HP\ = {v? : [0,T] x R rf measurable; < 00} and 

IMIfwi :== / / / / _ v \ 2H ~ 2 f(x ~ y)dydxdvdu. 

Jo Jo Js. d Jm d 

Remark 3.6. Using Fubini's theorem, we have the following alternative expression for calculating 
{ip, iP)hv'- for every ip,ip 6 2?((0, T) x R d ), we have 

(f,ip)nv = an I / \u-v\ 2H ~ 2 J 7 2<p(u,£)J : 2i>(u,QLi(d£)dvdu 

JO Jo JR d 

= a H I I I T 2 ^{u, i)\u - v\ 2H - 2 F 2 ^(v, £)dvdufjL(d£) 

Jg, d Jo Jo 



where 7ic(0, T) denotes the completion of 2?c(0, T) with respect to the inner-product (•, -)uc{o.t) defined 
similarly to (•, -)h(o,t)- particular, = J Rd H-T^yG, ?)IIh c (o rlM^O' This expression will not 

be used in the present paper. 

In this new context, the next theorem gives us a useful criterion for verifying that a function ip 
lies in TCP . To prove this theorem, we need the following lemma, generalizing Lemma 5.1, [llj . 

Lemma 3.7. For every A € Bb(R d ), there exists a sequence {g n ) n C £(M. d ) such that 

Proof: Let (</>„)„ C V(R d ) be such that (p n (£) — > 1a(£) uniformly and supp 0„ C K, Vn, where 
_K" C M. d is a compact (we may take 4> n = Ia * r) n , where rj n (x) — n dr q(nx) and 7/ £ 2?(R d ), with 
Jrii |*7(sc)]cfec = 1). Using the dominated convergence theorem and the fact that \i is locally finite, we 
get f Rd |^(0 - 1a(€)| 2 M^) - 0. 

Let V„ e S(R d ) be such that J" 2 ^n = 0n- Then 

Wn{i) - U(0lV(de) - 0. (25) 

Note that J Rd ^^(OlVC^O = / R d l&iCOlVf^) < 00 • By Lemma [CTTI ( Appendix C), it follows that 
VVi € 7 3 (IR d ). Since 7 ? (R <i ) is the completion of £ (R ) with respect to || • ||p(Rd), there exists a sequence 
(s-n)n C £(R d ) such that 

- ffnllpm,) = / Wn(0 " ^W£)| 2 /« - 0. (26) 

Jm d 

The conclusion follows from and (f!?6")) . ■ 
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Theorem 3.8. Let tp : [0, T] x M. d — > R be a function which satisfies the following conditions: 

(i) tp(-,x) € L 2 (0, T) for every x € R rf ; 

(ii) (t,x) h- > J-'o^TfiT, •) is measurable; 

^ T ^^^{T,x)f{x ~ y)T ^f{T,y)\T\'^ H ^ dydxdT < oo. 

T/ien (pel-CP. 



Proof: The proof follows the same lines as the proof of Theorem 12.21 The details are quite different 
though. Let 

A = : [0, T] x R d -> R; ^(-, x) £ L 2 {0, T) Mx, (r, x) h-> F >t<p(t, x) is measurable, and 
IMIi:=ci/ / / Fo,tv(t, x)f{x - r/)f ,rp(r, V )\t\- {2H -^ dydxdr < oo} 

JR JR d JR d 

A = {^:[0,r]xM^M; IMIi := ca C [ [ l"l 1/2 (u H ~ 1/2 ip(u, x))(s) ■ f(x - y) 

JO JR d JR d 

xI*: 1/2 (u H - 1/2 ip(u,y))(s) ■ s-^-^dydxds < oo} 
where c\ — cxhch and c 2 = {c* H T(H — 1/2)} 2 . The fact that 

Ac A and \\<p\\ K = \\<p\\a, Wtp £ A. (27) 
follows as in the proof of Theorem 12.21 let tp £ A is arbitrary; using the fact K* H is an isometry from 



W(0,T) to L 2 (0,T), we get (^(-, x), tp(-, y)) 2 H{%T) = (K* H tp(-, x), K* H tp(-, y))\ 2 ^ T) for all x,y £ R d . 
Multiplying by f(x — y), integrating with respect to dxdy and using Fubini's theorem, we get \\tp\\x = 
\\tp\\ A < oo. 

Next we prove that 

£ is dense in A with respect to || • ||a. (28) 
The proof of the theorem will follow from (|2"7]) and (|28p, as in the case of Theorem 12.21 

To prove ([28]). let tp £ A and e > be arbitrary. Let Xh(s) = s~( 2H ~^ds and 
a(s,x) — 1^1 l ^ 2 (u H ~^l 2 tp(u, x))(s). First, we claim that there exists g £ £ such that 

h ■■= I \ \ [a(s,x) - g(s,x)]f(x- y)[a(s,y) - g(s,y)]dydx\ H (ds) < e. (29) 
Jo Jm d JR d 

To see this, note that 

\\ l p\\a = c 2 / / a(s,x)f(x - y)a(s,y)dydxX H {ds) = c 2 / / |.F 2 a(s, £)|/x(d£)Afr(ds) < °°: 

JO JR d JR d JO JR d 

i.e. the map (s, £) i— * ^" 2 a(s,^) belongs to L 2 ((0, T) x R d ,<iA# x d/x) . Hence, there exists a simple 

function h(s,£) = E^i^T^-^.)^) 1 ^^) on (°> T ) x M ^ with ft e K > < 7j < ^ < T and B i c M<i 
Borel sets, such that 

T 

2 , 



\T 2 a{s,i) - h{s,i)\ 2 ^)\ H {ds) < e/4. (30) 

JR d 

By Lemma |3"771 for every j = 1, . . . , m, there exists gj £ £(R d ) such that 

/ \lB j (0-^g j (OHdO<^/(iD h ), (31) 

JR d 

where we chose D h = mY^ =1 P 2 \H{[lj,5j)). Define fl>(s,ir) = Y, 7 jLi 0j 1 h J ,s j )( s )9j{x). Clearly g £ £. 

Using ([317)) and (|3"Tj) , we get I\ — J Q T J Rd |.F 2 a(s, £) — ^gis, £)| 2 /z(<i£)Aff(ci,s) < e, which concludes the 
proof of f29]) . 
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We claim now that there exists a function I £ £ such that 

h -= / / [g(s,x) - b(s,x)]f(x - y)[g{s,y) - b(s,y)]dydxX H (ds) < e, (32) 



o 



where b(s,x) — 1^1 1 ^ 2 (u H ~ 1 / 2 l(u,x))(s). To see this, suppose that g — J2k=i b khc k ,d k ) i s )^A k (%) for 
some < Ck < dk < T and Ak C M. d Borel sets. By relation (8.1) of [T2), there exists an elementary 
function Ik £ £(0,T) such that 

'' T [l [ck , dk) (s) I%: 1/2 {u H -^l k {u)){ S )Y\ H {ds) < e/C g , (33) 



where we chose C g := || YZ=i kfc 1 ^JIp(R<«)- Let K s > x ) = Z)k=i b kh(t)lA k (x) £ £ and note that 

fc - =1 JO JR d JR d 

[l[ Cj , dj )(s) - lT: 1 ' 2 {u H - 1 ' 2 l J {u)){ S )]l Aj {y)dydx\ H (ds) 
k,j=l \ 9 9/ 

(we used (|3"3")l and the fact that a6 < 2(a 2 + b 2 )). The proof of (I32p is complete. 
Finally, we claim that from (|2U)) and (|3"2")1 , we get 

7:= / / [a(s,x) - b(s,x)]f(x - y)[a{s,y) - b(s,y)}dydx\ H (ds) < 4e, 

i.e. ||y — < 4ec2, which will conclude the proof of (|2"8"|) . To see this, note that J = 2fc=i -^fe> where 

h '= - I \ I [a(a,x)-g(s,x)]f(x-y)[b(s,y)-g(s,y)]dydx\ H (ds) 

JO JR d jR d 

h ■■= - I / / [a(s,y)-g(s,y)]f(x-y)[b(s,x)-g(s,x)]dydx\ H (ds). 

Jo JR d JR d 

The fact that I/3I < e and \1a\ < e follows from the Cauchy-Schwartz inequality in P(M. d ), (j2"5|) and 

(El. ■ 



It is again possible to describe the space TCP using the transfer operator. Define the transfer 
operator on £ by the same formula ([5]). Note that in this case we have 

( K H 1 [o,t]xA,KHli 0}S ] xB )-p = (^J K H (t, u)K H (s,u)duj (l A , l B )v(R d ) 

= R,H(t,s)(lA,lB)v(A d ) = (l[0,t]xA, l[0,s]xs)wP, 

i.e. K* H is an isometry between (£, (•, -)nv) and P ■ Since TCP is the completion of £ with respect to 
(•, -)rcPi this isometry can be extended to TCP . We denote this extension by K^ v . 

Lemma 3.9. K^ v : TCP — > P is surjective. 

Proof: The proof is similar to the proof of Lemma [2~3l using the fact that l[o,t]xA € ^^CHP) for all 
t £ [0, T],A £ B b (R d ), and £ is dense in P with respect to || • \\ v . ■ 
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Remark 3.10. Using ([9]) and Lemma [3T91 we can formally say that 

HV = {if such that (s,x) h-» fl-^-V^/^jVa^-i/a^^^jj^ is in v y 

3.2. The Noise and the Stochastic Integral. In this subsection, we introduce the noise which is 
randomly perturbing the heat equation. This noise is assumed to be fractional in time and colored in 
space, with an arbitrary spatial covariance function /. It has been recently considered by other authors 
(see e.g. [14]) in the case when / is the Riesz kernel and the spatial dimension is d = 1. The general 
definition that we consider in this subsection seems to be new in the literature. 

Let B — {B(ip); ip £ £>((0, T) x R d )} be a zero-mean Gaussian process with covariance 

E(B(<p)B(j>)) = (^i) HP . 

Let H B be the Gaussian space of B, i.e. the closed linear span of {B(tp); <p g 2?((0, T) x in L 2 (fl). 
As in subsection 12. 2i we can define B t (A) = -B(l[ ,t] X A) as the L2(^)-limit of the Cauchy sequence 
{B(ip)} n , where ((/?„)„ C T>((0,T) x R d ) converges pointwise to l[o,t]xA- We extend this definition by 
linearity to all elements in £. A limiting argument shows that 

E(B(<p)B&)) = { v , il))up, V<p, 4>e£, 

i.e. ip i — ► B(ip) is an is isometry between (£, (■, -)nv) an d H B . Since HV is the completion of £ with 
respect to (■, -)h, this isometry can be extended to HV, giving us the stochastic integral with respect 
to B. We will use the notation 

B(tp) = [ [ <p(t,x)B(dt,dx). 

JO JR d 

Remark 3.11. Similarly to subsection 12. 2i the transfer operator K^^, can be used to explore the 
relationship between B(tp) and another stochastic integral. Using Lemma |3.9[ we define 

M(0) := B{{K^ v )-\^>)), 0EV. (34) 

Note that 

E(M(<j>)M(r,)) = {{K^ v )- l {(t>),{K^ v )- l {ri))HV - (4,v)v, 
i.e. M = {M(cj)); <fi E V} is a Gaussian noise which is white in time and has spatial covariance function 
/. This noise has been considered by Dalang in [2]. We use the following notation: 

M{4>)= I [ cf>(t,x)M(dt,dx), 4>eV. 

Note that M(<f>) is in fact Dalang's stochastic integral with respect to the noise M. 

Let H M be the Gaussian space of M, i.e. the closed linear span of {M((j>); <f> & V} in L 2 (f2). 
By (EH), it follows that H M = H B . The following dia gram summarizes these facts: 



HV 



B 




B(cp) = M(K$i V <p), yep € HV, i.e. 



■/if Jr* v{t,x)B{dt,dx) = J T J Md (K^ v cp){t,x)M(dt,dx), e HV. 



H B =H M 
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In particular, B(t,A) = J* f A Kn(t, s)M(ds,dy). This relationship will not be used in the present 
article. 

3.3. The solution of the Stochastic Heat Equation. We consider the stochastic heat equation 
driven by the noise B, written formally as: 

Ut-Au = B, in (0,T) x R d , u(0, •) = 0. (35) 

As in subsection 12.31 we let Git, x) be given by formula fP3")) , and g tx (s, y) = G(t — s,x — y). The next 
theorem is the fundamental result leading to the necessary and sufficient condition for the existence of 
a process-solution and a distribution-solution of (|35p . 

To state the theorem we need to introduce the following notations: 



If,tx{r,s) ■= 9tx{s,y)f(y- z)g tx (r,z)dydz 

Jf(u,v,y,z) ■= G(u,y — x)f{x — x')G(v, z — x')dx dx' . 



Theorem 3.12. Suppose that the spatial covariance function f satisfies: 

A f (2t- s-r)- {d - a ?V 2 <I f , tx (r,s) < B f (2t-r-s)-^ d ~ a f^ 2 , (36) 

Vr e [0,i],Vse [0,t],VtG [0,T],VieR <l 
J f (u,v,y,z) < Cf(u + v)~ ( - d - a ?V 2 , (37) 
Vu£ [0,T],V»G [0,T],VyGM d ,VzeR d 

for some constants Af,Bf,Cf > and ctf < d. If 

H> d ~^, (38) 

then: (a) g tx £ \HV\ for every (t,x) € [0, T] x R d ; (b) rj * G £ KP for every rj £ V((0,T) x R d ). 
Moreover, 

WgtxWnv < oo V(t,af) G [0,T] x R d if and only if $M) holds. (39) 



Proof: (a) We will apply Theorem 13.81 to the function g tx . As we noted in the proof of Theorem 12.7 
9tx{-,y) £ L 2 (0,T) for every y £ R d , and the map (r, y) t— > !Fo,Tgtx(T,y) is measurable. We calculate 



htxllnv ■= a H c H 



Fo,T9tx(T,x)f(x - y)Fo,T9tx{r, y)\r\ {2H 1} dydxdr. 



For this, we write 



\gtxfnv = a H c H I \t\-( 2H -V I (£ e- iTS g tx (s,y)dsj f(y- z) ^ e iTr g tx (r,y)drj dy dzdr 
"//'// / |7-r< 2jr - x > / / e-^ s - r h fitx (r,s)drdsdr 



o jo 



t r t 

2H-2 , 



SB / / \s-r\ 2H -*I Ux {r,s)drds, (40) 
Jo Jo 
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where we used Lemma (Appendix) for the last equality and the definition of //^(r, s )- Using ([36]). 
we see that 

htAnv < a H B f f f \ S ~-r\ 2H - 2 {2t-r- S )- {d - a ^l 2 drds 
Jo Jo 

\\9t4nr > a H A f [ f \ S -r\ 2H - 2 {2t-r- S )^ d - a ^' 2 drds. 
Jo Jo 

Relation follows, since the integral above is finite if and only if 2H > (d — af)j2. In this case, we 
have ||5te|||«7>| = htxWnv < °o (since g tx > 0), and hence g tx S \HV\. 

(b) We will apply Theorem 13.81 to the function <p — rj * G. By Lemma T2.61 x) £ £2(0, T) 
for every x E M. d . We now calculate 



\<P\\nP := 



a H CH / / / ^o,Tf{T 7 x)f(x-y)J 7 ,Tf(r 7 x')\T\ (2H ^dx'dxdr. 
Jn JM d J~R d 



By (HSJ), we get 



M 2 nv = a H c H I \t\-v h -v r f T e —(s-r) 1 f r s r r v ( u+sMv+rjZ ) 

Jr Jo Jo JR d JR d Jo Jo 

Jf(u, v, y, z)dv du dz dy dr ds dr 

T r T r r rT-s r T-r 

2H-2 



= a H / / \s-r\ / / / / r)(u+s,y)r)(v + r,z) 

Jo Jo JR d JR d Jo Jo 

Jf(u, v, y, z)dv du dz dy dr ds, 

where we used Lemma [A. II (Appendix A) for the last equality and the definition of Jf(u, v, y, z). 

Using (|37|) and the fact that rj £ T>((0,T) x R d ) (and thus is bounded by a constant and its 
support is compact), we conclude that 

r T pT r r rT-s ?T- 

|2 s- „ n I I l„ „\2H-2 



MVnv < anC f / \ a -rY H - i / / / \n{u + B,y)ri{v + r,z)\ 

JO Jo JR d JR d Jo Jo 

(u + v)- {d - a i )/2 dy dz dv du dr ds 

rT rT rT-s rT-r 

< a H C f D v / \s-r\ 2H - 2 / (u + v)-^- a '^ 2 dv du dr ds. 

Jo Jo Jo Jo 

As in the proof of Theorem 12.71 (b), the last integral is finite if 2H > (d — a/)/2. 

■ 

The next theorem identifies identifies the constant a/ in the case of some particular covariance functions. 

Theorem 3.13. (i) If / = R a with < a < d, then and hold with atf = a. 
(ii) If f = B a with a > 0, then ffflj) and (37|) hold with a f = 0. 
(Hi) If f = G a with a > 0, then {5^1 and (57|) hold with a f = 0. 
(iv) If f — P a with a > 0, then and hold with af = -1 

Remark 3.14. a) If / = R a , condition (|38[) becomes H > max{(d — a)/4, 1/2}, which does not impose 
any restrictions on d. For any d > 1 arbitrary, it suffices to choose a such that max{d— 4iJ, 0} < a < d. 

b) If / = B a or / = G a , condition (j3"5)) becomes H > max{<i/4, 1/2}, which forces d < 4. 

c) If / = P a , condition (f3"5| becomes H > max{(d + l)/4, 1/2}, which forces d < 3. 
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Proof: We begin by examining condition 1136]) . To simplify the notation, we will omit the index tx in 
the writing of If,tx- Using the definitions of If and G, we obtain 



I f (r,s) 



1 



1 



(27T) 



/(V2(< - s)y' - y/2{t^r)z')e 



f(y-z)e 4 <*- s > 4 < t -') dydz 



dy dz 



= E[f{^2[T^s)Y-VW^r)Z)]=E[f{U)}. 



(41) 



Here we used the change of variables x — y — \/2{t — s)y', x — z = y/2(t — r)z' and we denoted by 
(Y, Z) = (Yi, . . . Yd, Zi, . . . , Zd) a random vector with independent N(0, 1) components, and 



U = y / 2(t-s)Y - y/2(t-r)Z. 



Note that U l = y/2(t - s)Y l - y /2(t-r)Z l , i = 1, . . . , d are i.i.d. 2V(0, 2(2t - s - r)) random variables. 
Then U 2 = Ui/y/2(2t - s-r), i = l,...,d are i.i.d. iV(0, 1) random variables and 

d d 

\u\ 2 =Y. u i= 2 ( 2t - s - r ) E ^ 2 = 2 ( 2< - s - r ^ ( 42 ) 



i=l i=l 

where is a xl random variable. 

We are now treating separately the four cases: 

(i) In the case of the Riesz kernel, f(x) = R a (x) = 7a, dM - '-' 4- "' 1 an d < a < d. Using (|^T|) 
and the integral If(r,s) becomes 

I Ra (r, a) = la , d E\U\- { - d -^ = la , d [m - s - r)]-^/ 2 E\Wd\- (d - a/2) 

:= C aA {2t - s - r)-( d - Q >/ 2 , (43) 

where C a , d = -/ aid 2~ {d ~ a)/2 E\ W d \~ {d ~ a)/2 ■ This proves that condition (H} is satisfied with a f = a. 

(ii) In the case of the Bessel kernel, f(x) — B a (x) — j' a f °° w lya ~ d ^ 2 ~ 1 er w e~' {X \~ l^dw and 
a > 0. Using (|4Tj) and (|42|) . the integral If(r,s) becomes 

/>oo 

Jo 

2f -r-s. 



exp 



2«i 



-W„ 



dw 



2t-r-s 



-d/2 



dw 



where we used the fact that E[e~ cWd ] = (1 + 2c) _d/2 for any c> 0. Note that 



2* - r - s 



w 



d/2 



< 1 



2t - r - s 



w 



d/2 



<c d 



1 + 



2t - r - s 



w 



d/2' 



where C d = 2 d / 2 ~ 1 . Hence 



„d/2 



— (2t-r- s)- d/2 < - 

2C d y ' ~ C d w d / 2 + (2f - s - r) d / 2 



< 1 



2t - r - s 



-d/2 



<w d l 2 {2t-s-r)- d l 2 (44) 



where for the first inequality we used the fact that a/ (a + x) > l/{2x) if x is small enough and a > 0. 
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We conclude that 

I Ba (r, s) < i a w^- d ^ 2 - 1 e~ w w d ' 2 {2t - r - s)- d ' 2 dw < i a T (|j (2t - r - s)~ d 2 



lB a (r, s) > 



( a -d)/2-i e - w ^ 2t _ r _ s yd/2 dw > _7a_ 

2C d 



•d/2 



2Cd Jo 

i.e. condition is satisfied with a/ = 0. 

(iii) In the case of the heat kernel, f(x) = G a (x) = 7" I ( Aa ) and a > 0. Using (141 [) and 

(|32"]l. the integral //(r, s) becomes: 

-<l/2 



Using (|44]) . we obtain that 



exp 



2t- s-r. 



2a 



7^1 + 



2t - r - s 



7^ 



(2t - r - s)- d / 2 < 7 Gq (r, s) < T « d a d/2 (2f - r - s)- d/2 , 



i.e. condition (|36|) is satisfied with a/ = 0. 

(iv) In the case of the Poisson kernel, f(x) = P a {x) = "f'^' d (\x\ 2 + a 2 )~( d+1 ^ 2 and a > 0. Using 
(|4T|) and (|42|) . the integral J/(r, s) becomes: 

(r, *) = lZ d E {(\U\ 2 + « 2 )- (d+1)/2 ] - lZ d E \2(2t - r - s)W d + a 2 f (d+1)/2 . 
Using the fact that 

A d [(2t -r- s)W d }- (d+1)/2 < [2(2t - r - s)W d + a 2 ]- {d+1)/2 < B d [{2t -r- s)W d }- (d+1)/2 
for some constants A d , B d > 0, we conclude that 

A d E\W d \- {d+1 ^ 2 {2t - r - s )-( d +D/ 2 < I Pa ( r , s) < B d E\W d \- {d+1 ^ 2 {2t - r - s)-^' 2 , 
i.e. condition (155]) is satisfied with ctf = — 1. 



We continue by examining condition f^Tf ). Using the definitions of J/ and G, we obtain that 

1 



J f (u,v,y,z) 



(4:TT) d {uv) d / 2 

1 



f(x — x')e 4 " 



4 " dxdx' 



(2ir) d 



f(y - z + V2(-</ua + yW))e 2 2 dada' 

= £[/(y-z + v^(v^y + v^JZ))] = £[/(y-z + t/)]. (45) 

Here we used the change of variables x — y = \[2ua and x' — z = \/2va', we denoted by (Y, Z) = 
(Yi, . . . Y d , Z\, . . . , Z d ) a random vector with independent N(0, 1) components, and 

U = V2uY - V2vZ. 

Note that Ui = \[2uYi — \[2vZi, i — 1, . . . , d are i.i.d. N(0, 2{u + v)) random variables. Then Vi = 
Uij \j2{u + v), i = 1, . . . , d are i.i.d N(0, 1) random variables and 

d d d 

\y-z + U\ 2 = J2(Vi ~ z i + ^) 2 = 2(u + v) + U,) 2 = 2(u + v) ^ if, (46) 
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where m — (j/j — Zi)j y/2(u + v) and Tj = fa + Vj, is N(/J,i, l)-distributed. It is known that (see e.g. p. 
132,0) 



(47) 



where Wd-i and S are independent random variables with distributions Xd-ii respectively N(y X^iLi Mf > !)• 
We are now treating separately the four cases: 

(i) In the case of the Riesz kernel, f(x) — r y a .d\x\^^ d ^ a ^ ■ Using (|43)l . (|4^|) . and (|47p . the integral 
J/(u, u, y, z) becomes: 



J Ra (u,v,y,z) = la , d E\y-z + U\- id - a ^ la .d[2(u + v)}- id - a ^ 2 E 



!=1 



-(d-a)/2 



= la , d [2{u + v)]-( d -^E \W d -i + S 2 \- (d - a)/2 
< D^u + v)-^/ 2 , 

where D a ,d = l a .,d2~ {d ~ a)/2 E\Wd-i\~^ a)/2 , i.e. condition (07D is satisfied with a f = a. 

(ii) In the case of the Bessel kernel, f(x) = B a (x) = 7^ J °° u;( Q ~ d )/ 2 ~ 1 e~" , e~' ;r ' 2 / 4 "'du; and 
a > 0. Using (gSJ), (gU), and (g7J) , the integral If(r,s) becomes 

J Ba {u,v,y,z) = i a / w^-^' 2 - l e- w E[e-\ y - z+u \ 2 /^]dw 



exp 



7^ / w^/^e^E 
Jo 

poo 



U + V 

2w 

U + V 

2w 



g?u> 



(W d _! + S 2 ) 



^[ e -^-i+S 2 )] < (i + 2c )- d/2 , Vc> 0. 



Note that 

This follows by the independence of W d -\ and S, the fact that E(e~ cWd - 1 ) = (1 + 2c)"( d ~ 1 )/ 2 , and 
£ (V^) = exp 1-%^ • i±|4 < (1 + 2c)" 1 / 9 , Vc > 



(48) 



(recall that 5 has N(y X^iLi M 2 > 1) distribution). Therefore 

J B >,W) < 7^ °^ (a - rf)/2 - 1 e- w (l + ^) 



-d/2 



dvdu 



= l' a w a/2 - 1 e- w {w + u + vy d ' 2 dw 
Jo 

< 7 ^r(a/2)( U + U )- d / 2 , 
i.e. condition (|57|) is satisfied with aj = 0. 
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(iii) In the case of the heat kernel, f(x) = G a {x) = i^ d e^ x]2 1 {Aa) and a > 0. Using g5]), gBJ), 
and (|4"T)) , the integral Jf(u, v, y, z) becomes 



exp 



u + v 
2a 



{W d -i + S 2 



< l" d 1 + 



u + f 



-d/2 



<7^> d/2 (- + «)- d/2 



where we used (l48l) for the first inequality. This proves that condition (|3T[) is satisfied with a/ = 0. 

(iv) In the case of the Poisson kernel, f(x) — P a {x) — ^ d {\x\ 2 + a 2 ) -( d + 1 )/ 2 and a > 0. Using 
(|45|) . (|46|) . and (|47|) . the integral Jf(u,v,y, z) becomes 

-(d+l)/2 



Jp>, W ) = 7 " d £ (\y- z + U\ 2 + a 2 )-( d+1 » 2 



2(u + v)J2T, 



2 , 2 



= 7^ |2(« + v)(W d ^ + S 2 ) + a 2 r {d+1)/2 < iMu + v)]-^y 2 E\W d ^\-^/ 2 



i.e. condition (|3~T)) is satisfied with 



Under the conditions of Theorem 13. 12|) . B(g tx ) and B(rj * G) are well-defined for every (t,x), 
respectively for every n £ 2?((0, T) x R d ), and we can introduce the following definition: 

Definition 3.15. a) The process {u(t,x);t € [0,T],x € R d } defined by 



u(t, x) := 5(5^) 



G(i — s,x — y)B(ds, dy) 



(49) 



/o Jm. d 

is called the process solution of the stochastic heat equation (I35p . with vanishing initial conditions, 
b) The process {u(rj); rj £ X>((0, T) x R d )} defined by 

u(r)) := B(rj * G) = [ [ ([ [ r](t + s,x + y)G(s,y)dyds] B(dt,dx) 

JO JR d \Jq JR d J 

is called the distribution- valued solution of the stochastic heat equation (|35|) . with vanishing initial 
conditions. 



Lemma 3.16. The process {u{t, x); t£ [0, T], x € R rf } is L 2 (VL)- continue 



Proof: The proof is identical to the proof of Lemma 12.101 based on the continuity of the function 
G(t,x) with respect to each of its arguments, and the fact that HfcH^-p < 00, which is a consequence 
of Theorem (a). ■ 

The next theorem can be viewed as a counterpart of the result obtained by Maslowski and 
Nualart (see Example 3.5, [8] in the case m = 1, L\ = A, / = 0, $ = 1). 

Theorem 3.17. Suppose that the spatial covariance function f satisfies 136\) and Let {u(i]);r) £ 

T>((0,T) x R d )} be the distribution-valued solution of the stochastic heat equation 1135)) . 

In order that there exists a jointly measurable and locally mean-square bounded process X — 
{X(t, x);te [0, T],x€ R rf } such that 

u(v)=[ [ X(t 7 x)n{t,x)dxdt Vry G P((0,T) x R d ) a.s. 
Jo Jwi d 
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it is necessary and sufficient that H38\) holds. In this case, X is a modification of the process u = 
{u(t,x);t£ [0,T},xeR d } defined by gPp. 

Proof: The proof is omitted since it is identical to the proof of Theorem l2.11[ using || • \\n-p instead of 
|| ■ ||ft, and relation (139)) instead of (fT5|) . ■ 



Appendix A. An Auxiliary Lemma 

The following result is the analogue of Lemma 1, p. 116, [15j . for functions on bounded domains. 
It plays a crucial role in the present paper. For our purposes, it is stated only for d — 1, but it can be 
easily generalized to d > 2. 

Lemma A.l. Let < a < 1 be arbitrary, (a) For every tp G L 2 (a,6) 7 we have 

b , 

|t|-( 1 -°V(*)* = Qc / \T\- a TaMT)dT 



where q a = (2 1 - a Tt 1 / 2 )~ 1 T(a/2)/T((l -a)/2). (b) For every ip, ip e L 2 (a,b), we have 
/ / <p{u)\u - v\- {1 - a) ^(v)dvdu = q a \T\- a T a , w {T)T a , b ^{T)dT 

Ja Ja Jr 

Remark A. 2. Note that q a — 1/7q,i where ^y a ,d is the constant defined in Example 13.1 



Proof: (a) We use the fact that 

f e -ir < e -^H 2 dT = 5 -i/2 e -\t\ 2 /(^S) i V <S>0. (50) 
Jr 

Using the definition of J- a ^(p, Fubini's theorem and (|50p . we have 

{ >„, 6 v(r)dT =1(1 e-^e-^dr) <p{t)dt = ^ 1/2 / '/(4»*> <p(t)dt. 

JR Ja \Jr J Ja 

Multiply by <5 Q / 2-1 and integrate with respect to 5 > 0. Using Fubini's theorem, we get 

/ OT 6a ' 2 ' le '^ T?d5 ) ?*Mr)dT = J (f°° s-^-oV^e-^/^dS^ <p(t)dt. 

Using the change of variable 1/5 = u for the inner integral on the right hand side, and the definition of 
the Gamma function for evaluating both inner integrals, we get the conclusion, 
(b) Note that for every u £ [a, 6], 

/ \u-v\-^- a ^(v)dv= f \w\-^- a ^(u-w)dw = q a [ |r|- Q jr u _ M _ (^)~(r)dT, 

Ja J u-b JR 

where we used the result in (a) for the last equality. Now, 

Fu-b,u-a{^u) \r) = f e- Tt il){u - v)dv = [ e- iT ^- w ^(w)dw = e-" u T a ^{r). 

J u—b J a 
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Using Fubini's theorem 



b />b r-b 

cp{u) / \u-v\~ ( ' 1 -^ip(v)dvdu = q a / if{u) I \T\- a e- iTU ^~^)drdu 



q a / / <p(u)e- lTU du \\T\- a T a ^(T)dr = q a / \T\- a T aib ip(r)^a,bHr)dT 



Appendix B. Proof of 
Let r) n (t,x) = A n (t- t )ijj n (x - x ) := a n (t)(3 n (x). Then 

T / /-O 







T , T {rin*G){T,x) = / e / / a n (t- s)/3 n (x-y)G(s,y)dyds ) dt 



Since supp a„ C (to, to + T/n), we obtain that 

/-T-s f Ft ,t +T/nUn{T) it S < t 

e- ir ^ s )a„(t - s)df = / e- l ™a„( U )d U - ^ ^_ MD+T/n a n (r) if t < -s < t + T/n 
J ~* [ if — s >t + T/n 

and hence 

^o,t(Vu * G)(t,x) = ^ . to+ T/ n anM / (3 n (x - y) e~ lTS G(s,y)dsdy+ 

JR d J-t 

e~' lTS G{s 1 y)J : ^ sMt+T / n a n {T)dsdy := A n (r,£c) + S„(r,x). 

-t -T/n 

Note that lim„ B n (r, x) = 0. Whereas for j4„(t, a;), we have 

nto+T/n pT/n 

Ft ,t +T/na n (T) = / e- lTt A„(i - to)* = e- lrt0 / e~ lTU \ n (u)du 

J t JO 

= e- lTta T .T\{T/n) -> e" irto , as n -> oo 

and 

&n(x-y) \ er lTS G{s 1 y)dsdy = \ ip n (x - y - x )J 7 - to ,oG(T,y)dy 



ip n {z - x )T-ta,oG(T, x - z)dz 

— > J--t ,oG(T, x — xq), as n — > oo 
by Lebesgue differentiation theorem (see Exercise 2, Chapter 7, [H]). Therefore 
lim A n (T,x) = e~" to T-t fiG{T,x - x ) = T^^gtoxo {r, x). 
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Appendix C. A Result about the Space V(R d ) 

The next result is embedded in Theorem 3, [2]. We have used this result in the proof of Lemma 
13.71 We include its proof for the sake of completeness. 

Lemma C.l. If ip E S(R d ) and J Rd |-F 2 ¥>(£)lV(dO < °°; then V e 

Proof: Let 77 £ T>(M. d ) be such that r\ > and J Rd r](x)dx — 1. Define rj n (x) = n d n(nx) and ip n = 
(p*r] n £ S(R d ). We have ip n G \V(R d )\ C V(R d ), since 



Wn{x)\f{x- y)\(p n (y)\dydx = / f(z)(\ip n \ * \(p n \){z)dz < 00, 
by Leibnitz's formula (see p. 13, [2]). Note that 

11^-^1^)= / i^n(o-^(oiV(de)= / ^(o-i^i^oim^-^o, 



where we used the dominated convergence theorem, and the fact that ^2^(0 = J~2rj{£,/n) — > 1 and 
1-^2^(01 — 1 f° r au Tl - The conclusion follows. 
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